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Abstract 

We calculate thermal corrections to Renyi entropies for free massless fermions on a 
sphere. More specifically, we take a free fermion on M x and calculate the leading 
thermal correction to the Renyi entropies for a cap like region with opening angle 20. 

By expanding the density matrix in a Boltzmann sum, the problem of finding the Renyi 
entropies can be mapped to the problem of calculating a two point function on an n 
sheeted cover of the sphere. We follow previous work for conformal field theories to 
map the problem on the sphere to a conical region in Euclidean space. By using the 
method of images, we calculate the two point function and recover the Renyi entropies. 



1 Introduction 


Entanglement entropy has become of interest to various communities in physics from con¬ 
densed matter and quantum information [TJ |2] to black holes and quantum gravity |3J, 0] . The 
discovery of a way of calculating entanglement entropy holographically produced interest in 
the AdS/CFT community as well [5j. 

In this paper, we use the conventional definition of entanglement entropy. We assume 
that the Hilbert space factors nicely with respect to two complementary spatial regions, A 
and A. The reduced density matrix and Renyi entropies are then defined as 


p A = tr A p, (1) 

S n = 3 —— log tr(p J 4) n . ( 2 ) 

1 — n 

The factor of 1/(1 — n) in the the definition of the Renyi entropy is convenient for taking a 
n —y 1 limit and recovering the entanglement entropy: 


See = ~ tr [p A log(p A )] = hin S, 

n —>1 


( 3 ) 


It was argued in Ref. [6] that for gapped systems at small temperature, thermal correc¬ 
tions to the entanglement entropy are Boltzmann suppressed. Further evidence in d — 1 + 1 
can be found in Refs. mMnnun]. For general conformal held theories with temperature 
1//3 on a circle of perimeter L , the coefficient of the Boltzman factor was calculated [12] : 
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SS n = S n (T) - S n ( 0) = 

1—71 

$See = See{T ) — See( 0) = 2gA 


n 


1 sin 2A (f) 

2A-1 0 ;„2A 


n 


Sill 


e -2nPA/L + O ( e -2K0A/L^ ( 4 ) 


7T^ f ITi 

1 — — cot ( — 


e -2^A/L + o(e -2^A/L )) 


( 5 ) 


where g is the degeneracy of the first excited state, A is the smallest scaling dimension 
among the operators, and i is the interval length. (In order for these formulae to hold, the 
conformal held theory on the circle has to have a unique ground state and a mass gap.) 

For higher dimensional conformal held theories on S 1 x § d_1 , an analogous thermal cor¬ 
rection to entanglement entropy is also known D3S- The result for the entanglement entropy 
for a cap on a sphere with polar angle 9 and radius R is given in general by the following 
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integral: 


SSee = gAI d (9)e-^ R + o (e~^ R ) , 


W) = 2vr 


Vol(s<M) f e dff c An : m w , in «y). 


( 6 ) 

(7) 


Vol^" 1 ) J o sin(0) 

The derivation of this result relies on a conformal transformation from the sphere to hyper¬ 
bolic space. Unfortunately, there can be subtleties associated with boundary terms when 
this transformation is invoked. In the case of conformally coupled scalars ns m, these 
boundary terms replace I d in the result above with I d _ 2 . The issue is that the conformal 
coupling requires a Gibbons-Hawking like term on the boundary. The natural constant 6 
boundary is different from the boundary of hyperbolic space, and this difference contami¬ 
nates the entanglement entropy. In this paper we study free fermions which have no such 
Gibbons-Hawking like term and consequently no subtleties associated with the boundary. 
Thus we expect and indeed find that the result ([b]) holds for massless free fermions. As an 
added benefit, we also compute thermal corrections to Renyi entropies for fermions. 

This paper is organized as follows. First we briefly review the mapping used in Ref. 
[133 that maps from the multi-sheeted cover of the sphere to a wedge in flat space. We 
then calculate the two point function using the method of images, from which we can read 
off the thermal corrections to Renyi entropies. Finally, we compare these calculations of 
entanglement and Renyi entropies with numerical results for fermions in d — 2 + 1 and 
d = 3 + 1. 


2 Renyi’s for a General CFT 


A main result from Ref. m was a general equation for the thermal correction to the Renyi 
entropy for a conformal field theory. We assume that when the conformal field theory is 
placed on a x M, there is a unique ground state |0) and a set of degenerate first excited 
states \ijji) with energy E^. We divide the into a spatial region A and complement A 
and consider instead of x M, an n-sheeted branched cover of this spacetime where the 
branching is over the region A. The result from Ref. [33] is 


ss n 


n ( {ifj i (z)jji(z')) n 

l ~ n i \(A(z)A(z'))i 



e ^ + o(e 


( 8 ) 
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where 'ip t is an operator that creates one of the first excited states. The point z is in the far 
Euclidean future and the point z! in the far Euclidean past. The subscript n indicates the 
two-point function is to be evaluated on this n-sheeted branched cover. (Note that the result 
Q could have been anticipated from a very similar result in 1+1 dimensional conformal held 
theories H90 

In general, it is not clear how to evaluate ('i/’i(z)'il>i(z')) n . However, if we restrict to the 
case where A is a cap on a sphere of opening angle 26, then we can take advantage of a 
conformal transformation that maps the n-sheeted branched cover of x K to C„ x 
where C n is an n-sheeted cover of the complex plane, branched over the negative real axis. 
It is convenient to make the transformation in a couple of steps, as was outlined in Ref. m- 
The Erst step takes the cap on to a ball in M d_1 (see the appendix of Ref. [15]): 


ds 2 = —dt 2 + dr 2 + r 2 dfl 2 

= fl 2 (—dr 2 + dd 2 + sin 2 ( 6 *)df 2 2 ) , 


where 


t ± r = tan 


t ± 9 


1 fr + 6\ ( t — 6 

\ l — - sec - sec - 

2 \ 2 V 2 


(9) 

( 10 ) 


( 11 ) 

( 12 ) 


and dfl 2 is the line element on 2 . If the cap has opening angle 26 o, then the ball has radius 

r o = tan($o/2). A further special conformal transformation maps the ball to a half space: 

= fl31 

J l-2i >-x+bV' [ ' 

ds 2 = dy^dy^d^ = (1 — 2 b ■ x + b 2 x 2 )^ 2 dx^dx u 8^, u . (14) 


We let x° and y° correspond to Euclidean times, and take b 1 = 1/ro to be the only non¬ 
vanishing value of the vector b. After further rescaling and rotations, the inserted operators 
can be placed at y' = ( 1 , 26 * 0 , 0 ) an d y = ( 1 , 0 , 0 ), where we are using polar coordinates (r, 9) 
on the C n . (For further details, see Ref. [14].) 

We will employ a method of images strategy for computing (^(y)^j(?/ / )) n on C n x R d ~ 2 . 
This strategy was already used successfully for the scalar in Refs. [HI [18]. The idea is 
to compute the two-point function using the method of images on the orbifold C/Z m for 
general m and then to obtain {^iiy^iiy')) n by analytic continuation, setting n = 1/m. As 
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the method of images relics on the fact that the underlying equations of motion are linear, 
we do not expect this method will be useful for interacting field theories. 

In the fermionic case, there are issues associated with nontrivial phases, signs and a choice 
of spin structure which we must address. One issue, which we now review, is that rotations 
act nontrivially on spinor wave functions. 

2.1 Rotation on Fermions 

For a Dirac fermion we know the effect of a rotation on the components of the spinor [TBj. 
A general Lorentz transformation in Euclidean signature, A, is given by 


1p(x) -> Al/2^(A 1 x), 


(15) 



(16) 


{l\ 7"} = 2 ST, 


(V) 


where parameterizes the rotations and Lorentz boosts. For the case of interest we are 
only interested in rotations in (0,1) plane, for which the only non-vanishing components are 
ujqi = — lo w = (j>- The matrix exponential can be done simply and is given by 


Ai/ 2 (0) = cos(0/2) + sin(^/ 2 ) 7 ° 7 1 . 


(18) 


If we then define 


y 2 = 7 ° + i 7 1 and 7 2 = 7° — *7 1 , 


(19) 


then equation (18) simplifies to 



( 20 ) 


A fact that we will rely on heavily moving forward is that 7 ° 7 2 and 7 ° 7 2 are projectors: 


(7V) 2 = 2( 7 V), (7VT = 2( 7 V), 

(tV)(tV)=0, (7°7 5 )(7V) = 0. 


( 21 ) 

( 22 ) 
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3 Analytic calculation of Renyi Entropies 


In flat space the fermion 2-point function is, up to normalization, 




7VO - z/% 

I y - y'\ d 

1 0 a d 1 

d — 2 |y — 2/'| d_2 


(23) 

(24) 


Following Ref. ra, the Green’s function on a wedge C/Z m x 2 can be calculated via 
the method of images. The Green’s function via the method of images is given by rotating 
one of the fermions by 2tt k/m where k indexes the wedges and m is the number of wedges. 
In going between adjacent wedges an extra factor of (-1) is added due to the spin structure 
(for example see Ref. [17]). The result is then 




7 < Vd At _ (-l) fc Ai /2 (27r/c/m) _ 

d — 2 ' ^ |j Z — / m ^|2 _|_ — y/^2j(d—2)/2 


(25) 


(Curiously, this expression only makes sense for m an odd integer. Nevertheless, we find that 
knowing the two-point function for odd integers is in general sufficient to make the analytic 
continuation to n = 1/m.) In the case of interest y = y' = 0, z' = e 2ie and z —>■ 1, the 
two-point function can be rewritten 




m—1 


— 7rik(m — 1) 7rifc(m — 1) — 

e ™ 7 "d z + e ™ 7 z dz 


k =0 


_ g2i(nk/m-\-0)\d —2 


7 

4(d-2) 


m—1 


V (P - 2) - id„) Y prr 


— 7rik(m — 1) 

6 m 


k=0 


^ _ g2i(7rk/m-\-Q)\d— 2 


+ 7 2 c.c. 


(26) 

(27) 


From this expression, we can deduce the following recursion relation for the two-point func¬ 
tion: 

Gf 1/ra , d+2) ( 20) = ((9, 2 + d(d - 2)bV + f) + 2t( 7 V - 7V)S.) (28) 

To obtain the Renyi entropy we make the replacement n = 1/m in the two point function 
and we use that f or f ree fermions on a sphere of radius R: 


SS n (9 ) = 


n 


1 — n 
n 

1 — n 


E 


(VW’i)] 


1 ) e -G~l)/3/(2R) + Q ( e -(d-l)0/(2R)j ^ 


tr(G^ id) (20)Gf ljd) (2d)- 1 - + 0 ( e -G-i)^/(2ip)_ 


(29) 

(30) 
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3.1 d=2 


In d — 1 + 1 we can choose gamma matrices (j° = a 3 and 7 1 


a 1 ) such that 


and 7 V = ( ° ° V 

1 ° 2 ) 

It is worth noting that it is convenient to have 7 ° 7 2 diagonal, but is not necessary. Then 
the 2-point function is given by 

C F 6XP( m ) | z eX P( m ) ) 

(1 /’ n,2) 2 ^ l 1 — exp (2i(kn/m + 6 1 )) 1 — exp (—2i(kir/m + 0)) / 

k—0 \ / (31) 

o ( z—is / g/ni jg . n\\ 

= 7 I 7 '—e esc(m 0 ) — j— e esc (mO) I . 

We can then calculate the Renyi entropies (and entanglement). 

sSn(t>) = hbte(Gf„, ! )(2«)Gf w) (29)- 1 - l)e-»« + o(<r«< 2s >) 

1 ^ (32) 

= —(sin(0) CSC (0/n) - n) + o(e~^ 2R) ) , 

SSee = 2 ( 1-0 cot( 0 ))e.- /3/(2R) + o(e^ /{2R) ) . (33) 



These agree with the known results for 2d CFTs [12] in general and for 2d fermions nun] 
in particular. 


3.2 d=4 


In d — 3 + 1 we choose gamma matrices 



In this case the 2-point function is given by 

ITT) — 

G'fi/ m 4 ) (20) = —(1 + 3 m 2 + (m 2 - 1 ) cos( 2 mO)) csc 3 (m 0 ) 7 ° - yV e ) . (34) 
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Repeating the calculation in d = 2 we get 


SS n (9) = r ^tr(Gf M) (20)Gf li4) (20)- 1 - l)e~ 3 ^ 2 ^ + o ( e -30/(2*)) 

= — 2 ((3 + n 2 — (n 2 — 1) cos(2 9/n)) esc 3 (9/n) sin 3 ($) — 4n 3 ) e ~ 3 /7(2-R) (35) 

+ o(e-WW) , 

SS ee {9) = lim -^tr(Gf n4) (20)Gf 14) (20)- 1 - i) e - 3 0/( 2 *) + 0 ( e - 3 /V( 2R )) 

n->i I — n K v 

= 2(5 + cos(2 9) - 60cot(0))e- 3/3/(2i?) + o{e~ zpl{ - 2R) ) . 

The second result correctly reproduces the entanglement entropy correction found for general 
conformal field theories on the sphere D3|- The result for 8S n (9) is new. 


3.3 d=3 


In odd dimensions we can choose the same as we would use in one smaller dimension. 
Namely, 

n l 2 0 \ n _ / 0 0 

= and 7°7~ = 

loo/ \ o 2 

Following previous work M [IS], we may try to convert the denominator of the Green’s 


function (27) to an integral in order to perform the sum over k. In the case of the scalar, 


the resulting expression can be analytically continued to all m and thus in particular to 
n = 1/m. However, in the case of the fermion, some extra phases appear to spoil the 
analytic continuation. We are able to extract thermal corrections to entanglement entropy 
from an n —>■ 1 limit of the integral successfully. Thermal corrections to Renyi entropies 
remain out of reach however. 

The first step in converting the sum to an integral is an integral representation of the 
cosecant used successfully in the analogous calculation for the scalars mm-- 


roo 0/'K+k/m—1 

/ dx — -= 7rcsc(9 + kn/m) . 

/q 1 + X 


From this integral representation, it directly follows that 


™} a —iTTk(m—l)/m 

E 


k =0 


sm(nk/m + 9) n 


"OO 0/71-— 1 ™ 1 

■4+7 E 

k =0 


k/m —iTck(m—l)/m 
X c- 


'0 


dx 


71 


x e/ir-l (3 + e -i*m x j 
1 + x 1 + e i ’ K / rn x 1 / m ' 


(37) 

(38) 

(39) 
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Using the representation (39) in the Green’s function (27) for d — 3, we obtain 


r F 
Wl/m, 


3) = g(l-^)7 Y- 


dx 


x 


6/ir—l p—inm(pim7r 


l + x 1 + e in / m x 1 /‘ 


e-"-" +x) 0 ^ 

+ 7 7 cc 


( 40 ) 


To get the entanglement entropy, we expand around m — 1: 

1 log(a: 


^(i/m,3) ~ (1 - ido) 


YY r°° 


dx 


9 

X*■ 


7r 


= 7V ie-* 0 csc 2 (0) 


l + i \a: 1 — x 

2e ie n 


(■m — 1) + 0(m — l) 2 + cc 


:{m — 1) + 0(m — l) 2 ) ) + cc . 


(41) 

(42) 


(1 + e ie ) 3 

The entanglement entropy correction is then constructed from a ratio of Green’s functions 

SSeeW = Jim Yl trtU,/™,3,(29)05, 3) ( 2 9)-‘ - l)<r Wii + o(e^' R ) (43) 

= 4tt csc(0) sin 4 (0/2)e- /3/i? + o(e~^ /R ) . (44) 


This result matches the general case derived in Ref. [11]. 

While this integral representation gives the correct thermal corrections to the entangle¬ 
ment entropy, it appears to fail for the Renyi entropies. We suspect a reason is that the inte¬ 
gral representation grows too quickly as a function of complex m to satisfy the assumptions 
of Carlson’s Theorem. In other words, there will not be a unique analytic continuation. We 
can break the integral up into two pieces, one from 0 < x < 1 and a second from 1 < x < oo, 
and then replace the two integrals with double sums: 


m— 1 


£ 


g— i7rk(m— l)/m 

sin(7rfc/m + 6) 


■£(-i ) p+ ” 

p,q=0 


m 


g— imitg— j7r(</-|-l)/m 


+ 


e -i-n(q+l)/m 


7r(l + nip + q) — m6 7r(l + m + nip + q) — m6 

ginq/m ^—imn ^inq/m 


n(mp + q) + mO n (m + mp + q) + md 


In the case of the scalar, the phases in the numerator of this expression vanish, and the sum 
has better convergence properties. Here instead, for m = iy pure imaginary, the sum has 
the same kind of growth as sin(7rm), which vanishes for all integer 


3.4 Recursion relation for Entanglement entropy 

We would also like to show that our recursion relation ( |28| is compatible with the recursion 
relation for the entanglement entropy found in Ref. HU- We start by Taylor expanding the 


1 See appendix A for an alternate integral representation of the sum. 
















two-point function and relating it to the entanglement entrop;^ 


G n , d (29) = G d (29) + 5G d (29)(n - 1) + 0(n - l) 2 , 

icsc d ~ 1 (6 ) 


G d {29 ) = 7°( 7 z e- ie - Ye ie ) 


2 d 


6S EE (6) = + o(e~“ E *). 


(45) 

(46) 

(47) 


We will proceed by induction and assume equation (|6]) in d dimensions. Equation (J7J) has 
the following recursion relation 

Vo\(S d ~ 2 ) sin d - 2 (0) 


UO) I d - 2 (e) 27r Vol(^-i) (d-l)(d-2)' 

Then using equation ([6]) and recalling that A — (d — l)/2, 

d — 1 


SG d (29 ) = 


2 

d — 1 


I d (9)G d (29) 

G d {29) (l d +2 + 2tt 


Vol(5 d ) sin d (0) 
Vol(S' d + 1 ) (d+l)dy ' 


Acting on both sides with the operator in equation (28) and simplifying yields 

SG i+2 ( 20) = ^tti/ J+ 2(9)G J+2 (29). 


(48) 

(49) 

(50) 

(51) 


We checked the entanglement entropy for both d — 1 + 1 and d — 2 + 1. Thus by induction 
the two recursion relations are in agreement for both even and odd dimensions. 


4 Numerical Check 

We are interested in numerically checking our results. As mentioned earlier a free fermion 
on a sphere does not suffer from the same boundary term ambiguities as the conformally 
coupled scalar 0- The numerics for a free fermion should then directly give the general con¬ 
formal held theory results. Using the convention T = ddy 0 the Hamiltonian and Lagrangian 
densities for a fermion in curved space are given by T 

C = (52) 

U = v^- // vI/(r,'7; ; Pl/. (53) 

(T a (x), ^ p (x')}y/ zr g = i8 a p8(x - x') (54) 

2 Here we are taking the case where 7 ° 7 Z is diagonal so that the inverse is particularly simple. 
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Where D\ is the covariant derivative on the manifold. This can be written explicitly in terms 
of the vierbcin (ej) and spin connection (cua/j)]^] We have defined the curved space gamma 
matrices and covariant spinor derivative 


Y = j'et , 


= 9/i + , 


(55) 

(56) 


where uj^ij is the spin connection. Using the torsion free Maurer-Cartan equation, de l Te- 7 A 
Uj l = 0, we can extract the spin connection, the nonvanishing elements of which are 


UJ 3 i = cos 


3 -1 

(°i) (n s[n ° k ) 

k=i +1 


(57) 


(58) 


The general Hamiltonian in d + 1 dimensions is then 

H (ll csc (^)) (&oi + ~Y~ cot ( 6 ^ vl; 

where i is a flat spatial index. 

We can remove the cotangents and the volume factor in the commutation relation with 
the following definition T = fllj=i csc V-Midj)) il>: 


{il3 a (x),^p(x')} = i5 a/3 8(x - x') , (59) 

d i —1 

u = n = $ 

i= 1 j =1 

To obtain a numerical result efficiently, we turn this Hamiltonian density into a d = 1 + 1 
Hamiltonian. To this end, we integrate over 6 t for i > 1. We then calculate the spectrum of 
O where the lowest energy, smallest eigenvalue, gives the lowest order thermal correction to 
the Renyi entropies. For general d the result is 

H d = [ d 6 ^ (7071 d 9l + ^ J 3 * * ( 61 ) 

Jo V 2sm{d 1 ) J 

as can be found in Ref. [19J. 

3 We use capital Roman letters J, J, K,... for flat space-time indices, lower case Greek for curved 

space-time indices, lower case Greek a, /3, 7 ,... for spinor indices, and lower case Roman *,j, k,.. . for curved 

spatial indices. 


7l<%i + 


sin(0i 




(60) 
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We discretize these Hamiltonians, turning the integral into a sum and the derivative into 
a finite difference. We then numerically calculate the entanglement and Renyi entropies in 
the same way as previous papers pa mug. We find agreement with our analytical results 
for the entanglement entropy in both d = 2 + 1 and d — 3 + 1 (see Figure [Tj) , and for the 
Renyi entropy in d — 3 + 1 (see Figure [2j) . 

«s ee ^/ rt 5S EE ^/ RT 




Figure 1: 5 See in d — 2 + 1 (left), 3 + 1 (right) with 200 grid points. 

tfpoE^ .SPQ^ 011 




Figure 2: For d — 3 + 1, SS 2 (left) and 5S 4 (right) with 200 grid points. 


5 Discussion 

In this paper we extended to include massless free fermions the work in Ref. na, which 
considered thermal corrections to Renyi and entanglement entropies for the conformally 
coupled scalar. This extension allowed a direct and successful comparison with Ref. [ 13 ] - 
which provided general results for thermal corrections to entanglement entropy for conformal 
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field theories - without additional complications caused by boundary terms present for the 
conformally coupled scalar. We also were able to calculate thermal corrections to the Renyi 
entropies for the free fermion in even dimensions. We give the analytic result in d — 1 + 1 


(32) and d = 3 + 1 (35) along with a recursion relation (28) which allows for computations 
of all even dimensions. In odd dimensions, we were unable to find an analytic continuation 
that would allow us to calculate thermal corrections to Renyi entropies, but we were able 
to reproduce the thermal corrections to entanglement entropy. Amusingly, the situation is 
usually reversed, where one can compute Renyi entropies but the analytic continuation to 
entanglement entropy is not feasible. 

It is possible that the methods used in this paper and those in Ref. 0 could allow for 
corrections to be calculated for other free higher spin theories and possibly more generally 
for conformal held theories. 


Acknowledgments 

We would like to thank Jun Nian and Kristan Jensen for useful discussions. We thank 
Ricardo Vaz for collaboration during the early stages of this project. We also thank Martin 
Rocek for help in finding the sum in Appendix A. C. H. and M. S. were supported in part 
by the National Science Foundation under Grant No. PHY13-16617. 


A Alternate Formulation of d=3 Sum 


In our effort to find the Renyi entropies in odd dimensions we came across an alternate form 


of the sum in Equation (38). 


™_1 p —mk(m—l)/m 

V —- 

sin(7T k/m + 9) 


= e- | oottntf) + i - 2-2— sill [ 2 ^ - */2m)1 

sm(mp) f~^/ 7 sm(7r k/m) 


(62) 


This alternate representation is essentially a Fourier series on the shifted interval < 9 < 
Using the same integral form for the cosecant used for d = 3 we can rewrite this new sum 
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as an integral. 


m/2—1 

E 


k= 1/2 


sin[2 k(6 — n/2m)\ 
sin(7r/c/m) 



dx 


rgk/m 

x(l + x) 


sin[2 k(9 


— 7 t /2 m)]. 


(63) 


The integral can be evaluated for individual values of n — 1/m. This rewriting of the sum 
seems to have different issues with analytic continuation than those that plagued Equation 
( |39j ). It appears to reproduce correctly the n — 2 thermal correction to the Renyi entropy 
(see Figure [3]), but fails for the others and for the entanglement entropy. 


6S 2 e“ R1 



Figure 3: For d — 2 + 1, SS 2 with 200 grid points. 


Equation (62) gives the following results for the first couple of Renyi entropies 

n = 2 2 sin(0/2) 3 

n = 3 | | [2 + cos (y)] s in 2 (|) 

(Note that the n — 3 result is not reproduced by our numerics.) 
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